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Complete evolution of the strongly interacting matter formed in ultrarelativistic heavy-ion col¬ 
lisions is studied within a coupled Boltzmann and relativistic viscous hydrodynamics approach. 

For the initial nonequilibrium evolution phase, we employ a multiphase transport (AMPT) model 
that explicitly includes event-by-event fluctuations in the number and positions of the participating 
nucleons as well as of the produced partons with subsequent parton transport. The ensuing near¬ 
equilibrium evolution of quark-gluon and hadronic matter is modeled within the (2+l)-dimensional 
relativistic viscous hydrodynamics. We probe the role of parton dynamics in generating and main¬ 
taining the spatial anisotropy in the preequilibrium phase. Substantial spatial eccentricities e n are 
found to be generated in the event-by-event fluctuations in parton production from initial nucleon- 
nucleon collisions. For ultracentral heavy-ion collisions, the model is able to explain qualitatively 
the unexpected hierarchy of the harmonic flow coefficients Vn(pr) (n = 2 — 6) observed at energies 
currently available at the CERN Large Hadron Collider (LHC). We find that the results for v u (pt) 
are rather insensitive to the variation (within a range) of the time of switchover from AMPT parton 
transport to hydrodynamic evolution. The usual Grad and the recently proposed Chapman-Enskog- 
like (nonequilibrium) single-particle distribution functions are found to give very similar results for 
v n (n = 2 — 4). The model describes well both the BNL Relativistic Heavy Ion Collider and LHC 
data for v n (pT ) at various centralities, with a constant shear viscosity to entropy density ratio of 
0.08 and 0.12, respectively. The event-by-event distributions of 02,3 are in good agreement with the 
LHC data for midcentral collisions. The linear response relation v n = k n £n is found to be true for 
n = 2,3, except at large values of e n , where a larger value of k„ is required, suggesting a small 
admixture of positive nonlinear response even for n = 2, 3. 

PACS numbers: 25.75.Ld, 24.10.Nz, 12.38.Mh 


INTRODUCTION 

High-energy heavy-ion collision studies at the BNL 
Relativistic Heavy Ion Collider (RHIC) fil 0 and the 
CERN Large Hadron Collider (LHC) E-fy] have firmly 
established the formation of a strongly interacting quark- 
gluon plasma (QGP) close to local thermodynamic equi¬ 
librium. Evidence for this is provided by the hydrody- 
namical analyses of collective flow, which require an ex¬ 
tremely small shear viscosity to entropy density ratio p/s 
0 El- The precise extraction of p/s depends crucially 
on the adequate knowledge of the initial-state dynam¬ 
ics, the lack of which, at present, represents the largest 
uncertainty in the hydrodynamic modeling of heavy-ion 
collisions 

Smooth (nonfluctuating) initial energy-density distri¬ 
butions based on the Glauber or color-glass-condensate 
(CGC) models have been used successfully to describe 
the elliptic flow i; 2 , albeit with differing values of p/s 
{l3| . On the other hand, event-by-event fluctuations in 
the initial configuration of the system are solely respon¬ 
sible for the (rapidity-even) odd flow harmonics, V 3 , v$ 
0- The fluctuations are also important in explaining 
the double-peaked structure in the finaldihadron corre¬ 
lations across large rapidities SSElEEI. 


The Monte Carlo method has been invoked to simu¬ 
late event-by-event geometric fluctuations in the number 
and positions of the participant nucleons (MC-Glauber) 
@-12 1 or of intrinsic gluons (MC-CGC) [3- Both meth¬ 
ods give rise to “lumpy” initial conditions and are able 
to reproduce elliptic flow data at RHIC and LHC EE1- 
MC-CGC yields a relatively larger initial eccentricity and 
and hence requires larger viscous damping (p/s) to ex¬ 
plain the elliptic flow data . However, it underpredicts 
the triangular flow data 17]. Most of these model stud¬ 
ies “switch on” hydrodynamics at a specified early time, 
ignoring the preequilibrium dynamical evolution towards 
the hydrodynamic regime. Indeed, various ansatze and 
parametrizations have been adopted to describe the ini¬ 
tial distribution of the unknown kinetic and thermody¬ 
namic quantities in these models. 

In a more recent study, the preequilibrium dynam¬ 
ics was approximately accounted for by combining an 
impact-parameter-dependent saturation model with the 
classical Yang-Mills description of the flowing glasma 
fields 18J. The model considered the color charge fluctu¬ 
ations inside the nucleon as well and could describe re¬ 
markably well the transverse-momentum (pr)-dependent 
and pT-integrated flow harmonics v n (n = 1 — 5) data 
at LHC @, |J]. A fully dynamical simulation which in- 






2 


eluded solutions to anti-de Sitter/conformal field theory 
(AdS/CFT) for the preequilibrium stage and the stan¬ 
dard viscous hydrodynamics for the equilibrium stage 
was recently developed 0- Preequilibrium dynam¬ 
ics has also been accounted for in the ultrarelativistic 
quantum molecular dynamics (UrQMD) string dynamics 
model [§]. These studies underscore the importance of 
including the inherent fluctuations at both nucleonic and 
partonic levels as well as the initial-state preequilibrium 
dynamics. 

The aim of the present paper is to study various as¬ 
pects of the hydrodynamic flow in a model with an initial 
state that includes fluctuations at the nucleonic as well as 
subnucleonic levels and which moreover considers pree¬ 
quilibrium dynamics of partonic matter. For this purpose 
we have used a multiphase transport (AMPT) model 
[0 that consistently incorporates these fluctuations in 
the preequilibrium phase. For the subsequent near¬ 
equilibrium evolution, the standard (2+l)-dimensional 
[(2+l)D] viscous hydrodynamics code VISH2+1 fj] is 
used until the matter reaches freeze-out. 


THE MODEL 


The version of AMPT used in this paper has the initial 
conditions based on the HIJING 2.0 model [2lJ,[22j, which 
employs the Glauber model with Woods-Saxon nuclear 
distribution to determine the nucleon configuration in an 
event. While the soft (binary) nucleon-nucleon collisions 
lead to string excitations, the hard collisions produce un¬ 
correlated minijet partons. Fluctuations in the minijet 
parton multiplicity follow the Poisson distribution where 
the average multiplicity (determined by the minijet cross 
section) increases strongly with the center-of-mass en¬ 
ergy [23[. We have employed the string-melting version 
of AMPT 0] for the initial partonic state, which consis¬ 
tently incorporates these fluctuations and describes bet¬ 
ter the final collective behavior of hadrons. In this version 
the strings melt into their constituent quarks and anti¬ 
quarks, whose positions are determined randomly within 
the string, which is yet another source of fluctuations. 
The scatterings among these quarks and minijet partons 
are treated with Zhang’s parton cascade (ZPC) with a 
parton-parton elastic cross section of 1.5 mb. 

At the switchover time r sw , from the microscopic 
AMPT model to the macroscopic viscous VISH2+1 hy¬ 
drodynamic model [7j, the local energy-momentum den¬ 
sity, baryon number density, and the flow velocity in an 
event can be estimated from the positions and momenta 
of the formed partons at midrapidity: We use the space- 
time rapidity window [—1,1]. To remove the possible 
numerical instabilities in the ensuing hydrodynamic evo¬ 
lution, each pointlike parton in AMPT is smeared with 
a 2D Gaussian distribution in the transverse plane 


( a ) b = 0 fm 



FIG. 1: (Color online) Energy-density distribution in the 
AMPT model at the switchover time r sw = 0.4 fm/c in the 
transverse xy plane, with the smearing parameter a = 0.8 
fm. The results are for (a) central (6 = 0 fm) and (b) noncen¬ 
tral (6=5 fm) Pb+Pb collisions at a center-of-mass energy 
%/snn = 2.76 TeV. 


Energy density in the local rest frame then becomes 
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where ( Xi , yi) are the transverse coordinates and E[ is 
the energy of the ith parton in the local rest frame. The 
parameter N is introduced to account for the various 
uncertainties inherent in the present approach, such as 
those owing to our prescription to match the (3+l)D mi¬ 
croscopic description with the (2+l)D macroscopic one, 
at the selected switchover time r sw and the smearing pa¬ 
rameter <t, any other missing preequilibrium physics, lack 
of hadronic after-burner, etc. The value of N was ad¬ 
justed to describe the most central px spectra. N was 
found to be of the order of unity: 1.38 at RHIC and 1.24 
at LHC. The same value was used at all other centrali¬ 
ties. A similar normalization constant has been used in 
Refs. [ 1 , 0 , [ 3 , 0 ] • The Gaussian width is set at a = 0.8 
fm. Larger (smaller) values of a lead to smaller (larger) 
fluctuations and hence reduced (enhanced) magnitudes 
of the odd flow harmonics, v n . Final numerical results 
were found to be insensitive to the choice of the initial 
transverse flow velocity 0]; hence, we have set it equal 












3 


to zero. A detailed study of the effects of the initial flow 
velocity was performed in Ref. 24]. 

Figure [T] shows the energy-density profiles in the trans¬ 
verse plane, calculated in the AMPT model, for Pb+Pb 
collisions at yjS nn = 2.76 TeV, in a central (impact pa¬ 
rameter 6 = 0 fm) and a noncentral (6 = 5 frn) event. 
The snapshots are taken at the switchover time r sw = 0.4 
frn/c. As expected, the high-energy-density zone is of 
wider spatial extent in the central collision than in the 
noncentral collision. We further notice that these dis¬ 
tributions are azimuthally anisotropic and also contain 
several peaks which correspond to local maxima or “hot 
spots” [251 ] . which result in enhanced pressure gradients. 

The hydrodynamic evolution is governed by the con¬ 
servation equations for particle current, = 0, and 

the energy-momentum tensor, d^T^ 11 = 0, where 


= nu^ + 

= eu»u v - PA ^ + 7. (2) 




FIG. 2: (Color online) Transverse-momentum spectra of 7r + , 
K + , and p for two centrality ranges, 0 — 5% and 20 — 30%, in 
Pb+Pb collisions at y /snn = 2.76 TeV. The symbols repre¬ 
sent ALICE data 0] at midrapidity and the lines correspond 
to AMPT+Hydro calculations. 


Here P, n, e are, respectively, hydrodynamic pressure, 
number density, and energy density, and = g^ — 
rPw" is the projection operator on the three-space or¬ 
thogonal to the hydrodynamic four-velocity u M defined in 
the Landau frame: = eit M . The dissipative quan¬ 

tities, viz., the particle diffusion current and the shear 
pressure tensor, are denoted by rd 1 and , respectively. 
In this work, we set the net baryon number current, the 
bulk viscous pressure, and the initial shear viscous tensor 
to be zero. 

We have employed the s95p-PCE equation of state 
(EoS) [26|, which is obtained from fits to lattice data 
for crossover transition and matches a realistic hadron 
resonance gas model at low temperatures T, with partial 
chemical equilibrium (PCE) of the hadrons for tempera¬ 
tures below XpcE ~ 165 MeV. 

The hydrodynamic evolution is continued until each 
fluid cell reaches a decoupling temperature of Idee = 120 
MeV. The hadronic spectra are obtained at this temper¬ 
ature using the Cooper-Frye freeze-out prescription 271 

= <3) 

where Y is the rapidity, g is the degeneracy factor, 
is the particle four-momentum, dY,^ represents the el¬ 
ement of the 3D freeze-out hypersurface and f{x,p) is 
the nonequilibrium phase-space distribution function at 
freeze-out, which can be written as a small deviation from 
the equilibrium distribution function /o, i.e., / = fo + Sf. 
Unless otherwise mentioned, we use the viscous correc¬ 
tion form corresponding to Grad’s 14-moment approxi¬ 
mation f28| . 

6f= 2 (e + P)T* pap ^ (4) 


where corrections up to second order in momenta are 
present, and f 0 = l — rfo, with r = 1,—1,0 for Fermi, 
Bose, and Boltzmann gases, respectively. Resonances of 
masses up to 2.25 GeV are included in the calculations 
(so as to be consistent with the s95p-PCE EoS). Results 
presented here include resonance decays; this tends to 
reduce the anisotropic flow especially at low transverse 
momenta. 


RESULTS AND DISCUSSIONS 

Figure [2] shows the transverse-momentum spectra of 
pions, kaons, and protons in the VISH2+1 calculation 
with AMPT fluctuating initial conditions in the 0 — 5% 
and 20 — 30% central Pb+Pb collisions at yj snn = 2.76 
TeV in comparison with the ALICE data at midrapidity 
fiol . To connect the centrality c to the impact parameter 
6, we use the empirical relation c = 7 t6 2 /ct ,30], with the 
nucleus-nucleus total inelastic cross section a = 784 fm 2 
calculated from the Glauber model. The results are ob¬ 
tained with 200 events, with a shear viscosity to entropy 
density ratio rj/s = 0.12 in the VISH2+1 calculation. 
Here and in subsequent figures, unless stated otherwise, 
error bars for the relevant pt range are comparable to 
the symbol sizes. We find that the spectra for 7r + and 
K + from the full AMPT+Hydro simulations are in good 
overall agreement with the experimental data. However, 
the proton yield is overpredicted in the model possibly 
owing to the neglect of the hadron spectrum above 2.25 
GeV and final-state hadronic rescatterings, particularly 
pp annihilations. Incidentally, this proton excess is in 
line with the “proton anomaly” seen in the statistical 
hadronization model fits to the LHC data at this energy 

31|- 
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FIG. 3: (Color online) Time evolution of event-averaged 
eccentricity (En) and px dependence of root-mean-square 
anisotropic flow coefficients (a 2 ) 1 / 2 for charged hadrons in 
the AMPT+Hydro calculations for head-on (b = 0) Pb+Pb 
collisions at sj snn = 2.76 TeV. The vertical line at r = 0.4 
fm/c corresponds to the switchover time. Results are for the 
nucleon distributions in the colliding nuclei, which are az- 
imuthally asymmetric (a),(b) and symmetric (c),(d); see text 
for details. The symbols represent the CMS data [34} for 
u n {2} obtained from two-particle correlations in the 0-0.2% 
ultracentral Pb+Pb collisions at LHC. 


The nth harmonic participant eccentricity in the 
AMPT model can be estimated from the transverse po¬ 
sitions of the produced partons as 32] 


J(r n cos (ruj))) 2 + (r n sin (ruj))) 2 , 

Sn = -<r”>-' (5) 


Here (r, <f>) are the polar coordinates of the partons (or 
of the nucleons at r = 0 ) in the transverse plane, and 
the average (• • •) in a given event is obtained by using 
the smeared local energy density as the weight. For the 
calculation of e n during the lrydrodynamical evolution, 
(r, (j>) refer to the cell coordinates. Figure ©a) depicts 
the time evolution of event-averaged eccentricities in cen¬ 
tral (b = 0) collisions. For these collisions, eccentricities 
arise purely from event-to-event fluctuations of particle 
positions. The higher-harmonic eccentricities are found 
to drop faster with time for r > r sw = 0.4 fm/c. The 
minima correspond to the change in the sign of e n at 
late times. 

We recall that the fluctuations arise at both nucleonic 
and partonic levels. To isolate the contribution of par- 
tonic fluctuations, we consider the nucleon configuration 
obtained via Monte Carlo Glauber in one quadrant of the 
overlap zone and replicate it in all other quadrants, by 
imposing reflection symmetry with respect to the x and y 
axes, in each event. This results in a nucleon distribution 
which is approximately azimuthally symmetric in the ini¬ 
tial state. As seen in Fig. EKc), this causes eccentricities 


(e n ) to vanish at r = 0. However, the subsequent pro¬ 
duction of partons through (semi-)hard nucleon-nucleon 
collisions results in a sizable and rapid generation of (e n ). 
For even harmonics (not shown in the figure) the rise is 
higher than that for odd harmonics, owing to the initial 
azimuthal symmetry. 

The flow harmonics v n are generated by conversion of 
the initial spatial anisotropy into momentum anisotropy 
owing to unequal pressure gradients in the hydrodynamic 
evolution. To calculate the v n in the event-by-event fluc¬ 
tuating AMPT plus VISH2+1 hybrid model and compare 
it with the LHC data, we first estimate for each harmonic 
the event-plane angle 'F„ as in Emm 


4b 


1 

= — arctan 
n 


(sin {ruj))) 
(cos (ruj))) ’ 


( 6 ) 


where (f> is the azimuthal angle of the outgoing hadron 
momentum. The flow coefficients are then determined 
using 


v u (pt) = (cos [n(<j> - 4+)]} 

_ f d(f> cos [n((j> — 4 , „)] dN/(dY px dpx d(j>) 

J d(j> dN/ (dY pt dpT d(j>) 

(7) 

where dN/(dY px dpx d(j>) is the particle spectrum. As 
in Eq. the angular brackets in Eqs. © and © also 
refer to averaging over all the relevant particles in a given 
event. In the present (2+l)D (longitudinal boost invari¬ 
ant) viscous hydrodynamic simulation, the dependence 
on the rapidity Y is suppressed. We calculate the root- 
mean-square value of the flow harmonics, (v^Px)) 1 ^ 2 , 
over the entire event sample and compare it with the v n 
data at RHIC and LHC. We recall that for the v n {2} data 
obtained with the two-particle correlation method or the 
v n {EP} data obtained with the event-plane method, the 
root-mean-square (rms) y n is the most appropriate quan¬ 
tity to compare with [331 ]. 

Figures ©b) and[3][d) show the px dependence of the 
rms v n for charged hadrons at small px- For the az¬ 
imuthally symmetric nucleon distribution, £3 and €5 origi¬ 
nating purely from parton fluctuations translate into siz¬ 
able fractions (~ 50%) of the total ^3 and v 5 , respec¬ 
tively; see Fig. ©d). 

Figure[3]]b) also displays an intriguing feature observed 
in the CMS data for v n (px) in ultracentral (0 — 0.2%) 
Pb+Pb collisions at sJsnn = 2.76 TeV: Elliptic flow 
V 2 (pr) has the largest magnitude up to px — 1 GeV/c 
and thereafter it flattens out, whereas the higher har¬ 
monics continue to rise and become successively larger 
than V 2 with increasing px■ This is in contrast to what is 
observed at higher centralities; see Fig. 0 for instance. 
Our hybrid-model calculations are in qualitative agree¬ 
ment with the CMS data over the entire px range. Better 
understanding of ultracentral collisions is still a challenge 
for the theory 35|, 
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FIG. 4: (Color online) Root-mean-square anisotropic flow 
coefficients, (u^(pt)) 1 ^ 2 (for n = 2-6, top to bottom), for 
charged hadrons in the AMPT+Hydro calculations in 20 — 
30% central Pb+Pb collisions at yj snn = 2.76 TeV and 
r//s = 0.12. The results are (a) for two different switchover 
times, T sw = 0.4 fm/c (solid lines) and r sw = 0.8 fm/c (dashed 
lines); and (b) for two different forms of Sf in the freeze-out 
prescription, viz., Grad’s approximation [Eq. Q] (solid lines) 
and Chapman-Enskog expansion [Eq. ([8])] (dashed lines) with 
r sw = 0.4 fm/c. 



FIG. 5: (Color online) Transverse-momentum dependence of 
the rms anisotropic flow coefficients (u^(pt)) 1 ^ 2 of charged 
hadrons calculated at various centralities in Pb+Pb collisions 
at y/SNN = 2.76 TeV in the AMPT+Hydro hybrid approach 
(lines) with 77 /s = 0.12 as compared to the ATLAS data 
(symbols) based on the event-plane method. 


Figure HD a) shows a comparison between v u (pt) ob¬ 
tained using two different switchover times, r sw = 0.4 
fm/c (solid lines) and r sw = 0.8 fm/c (dashed lines), for 
20 — 30% central Pb+Pb collisions at y/S nn = 2.76 TeV 
at a fixed 77 /s = 0.12. We find that within a window of 
0.3 < r sw < 0.8 fm/c, a larger switchover time slightly 
reduces v n (j>T ), especially for lower harmonics. However, 
the difference is negligible, which suggests that the flow 
buildup at very early times is rather insensitive to the 
transport or hydrodynamic modeling. This also indicates 
that the parton dynamics in AMPT quickly brings the 
system to local equilibrium. The lack of sensitivity to 
the switchover time was also noted in Refs. H®. 

We now explore the uncertainty in the flow owing 
to a different choice of the viscous correction in the 
nonequilibrium distribution function. An alternate form 
to Grad’s approximation, Eq. d4]), is based on the 
Chapman-Enskog-like approach [37(, 


Sf 


5/o/q 

8 PT(u-p) 




( 8 ) 


which is obtained by iteratively solving the Boltzmann 
equation in the relaxation-time approximation. It may 
be noted that, though the two nonequilibrium distri¬ 
butions lead to distinct viscous corrections, the dissipa¬ 
tive evolution equations derived from them have identical 
forms and coefficients in the extreme relativistic limit. 
In Fig. |U[b), we compare v u (pt) based on Grad’s ap¬ 
proximation (solid lines) and Chapman-Enskog expan¬ 
sion (dashed lines) for the same initial conditions. The 
two forms of Sf result in essentially similar v n within 


the statistical errors. Note, however, that the Chapman- 
Enskog results tend to be slightly below Grad’s at small 
Pt and slightly above Grad’s at larger pr- This is at¬ 
tributable to the fact that the Chapman-Enskog method 
leads to essentially linear momentum dependence of Sf as 
opposed to the quadratic dependence in case of Grad’s 
approximation. This results in larger (smaller) viscous 
corrections in the Chapman-Enskog approach at smaller 
(larger) pt 37\. In the following, we use r sw = 0.4 fm/c 
for switchover time and Grad’s approximation for Sf in 
the freeze-out prescription. 

Figure [5] shows our results for the rms v u (pt), in 
comparison with the ATLAS data [I] obtained in the 
event-plane method, at various centralities. We find 
overall good agreement with the data for all harmonics 
(n = 2 — 6 ) and at all centralities. A single fixed value 
77 /s = 0.12 of the shear viscosity to entropy density ra¬ 
tio is able to achieve the required suppression (relative 
to perfect hydrodynamics) of the flow harmonics for all 
centralities. Apart from 77 /s, flow also depends on the 
smearing parameter <7 appearing in Eq. JT|). It controls 
the granularity and hot spots in the initial state [11]. A 
smaller a results in a hardened pt spectrum owing to 
larger pressure gradients [§|. Elliptic flow u 2 is less sen¬ 
sitive to variations in a , compared to the triangular flow 
V 3 , which is entirely driven by fluctuations. A smaller 
value of er tends to raise the value of U3, which can be 
compensated by choosing a higher value of 77 /s. However, 
this would lead to underprediction of u 2 . This suggests 
that the initial-state fluctuations alone are unlikely to re¬ 
produce all the observed w„, and the interplay between 
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FIG. 6: (Color online) (a) Centrality dependence of the rms 
values of the anisotropic flow coefficients (u^) 1 ^ 2 of charged 
hadrons in Pb+Pb collisions at yj snn = 2.76 TeV calcu¬ 
lated in the AMPT+Hydro (solid lines) with 77 /s = 0.12 as 
compared to v n {2} obtained from two-particle correlations 
within the complete AMPT calculations (dashed lines) and 
the ALICE data Q] (symbols), (b) Centrality dependence 
of the rms values of the eccentricities (e 2 ) 1 ^ 2 at the hydro- 
dynamic switchover time, (c) Centrality dependence of the 
ratio (u 2 ) 1,/2 /(e 2 ) 1/,2 in the AMPT+Hydro model. 


fluctuations and viscosity is essential. As several physical 
sources of fluctuations in both the nucleons and partons 
contribute to the AMPT initial conditions (as already de¬ 
scribed), a reasonably diffused initial state with er = 0.8 
fm is required to explain the v n (n = 2 — 6) ATLAS data. 

In Fig. 0a), p-r-integrated rms values of v n in AMPT 
plus VISH2+1 calculations are compared with ALICE 
data 0 for v n {2} obtained from two-particle correla¬ 
tions, as a function of centrality. Our calculation of 
the pT-integrated flow takes into account the Jacobian 
appearing in the transformation between rapidity and 
pseudorapidity 38]. It also takes into account pseudo¬ 


rapidity and pr cuts appropriate for the ALICE data. 
With the same constant ij/s = 0.12 at all centralities, 
the model shows a remarkable agreement with the data 
for all harmonics. As expected, V 2 exhibits a strong cen¬ 
trality dependence as it is driven mostly by the initial 
spatial anisotropy, whereas the odd harmonics are en¬ 
tirely induced by fluctuations. 

Figure El a) also shows the v n {2} obtained solely in 
the AMPT model without the inclusion of VISH2+1 hy- 



FIG. 7: (Color online) (a),(b) Correlations between harmonic 
flow coefficients and eccentricities at r = r sw for n = 2 — 3. 
The solid lines represent the fits, (c) Event-by-event distri¬ 
butions p(v n ) calculated in AMPT+Hydro with p/s = 0.12 
(solid lines) and the eccentricity distributions, p{s n ) scaled 
to the respective (v n ) from AMPT initial conditions (dashed 
lines). Errors on the calculated distributions are shown. Sym¬ 
bols represent the ATLAS data [4lJ. Results are for charged 
hadrons at pr > 0.5 GeV/c in 20 — 25% central Pb+Pb colli¬ 
sions at yj snn = 2.76 TeV. The analysis was performed with 
a sample of 3000 events. 


drodynamic evolution. This is slightly larger than that 
obtained in the AMPT+Hydro calculation. This can be 
understood by noting that the hybrid calculation employs 
a softer EoS with a crossover transition and a larger vis¬ 
cous damping, relative to the pure AMPT calculation 

BIH. 

Figure 0b) shows the rms value of the eccentricity at 
r sw = 0.4 fm/c, as a function of centrality; see Eq. 10. 
In the limit of ultracentral collisions, the eccentricities 
arise solely from fluctuations and are found to be nearly 
identical. As the impact parameter increases, they all 
increase: The ellipticity becomes more pronounced for 
midcentral collisions, and for the higher harmonics the 
fluctuations become more important for smaller overlap¬ 
ping geometry. However, as expected, £2 > £3 ~ £4 « £5 
as the shape of the overlap zone is predominantly elliptic. 

We now study the efficiency of conversion of the initial 
spatial anisotropy into the final momentum anisotropy 
as a result of hydrodynamic evolution. To that end, we 
present in Fig. 0c) the ratio k n of rms values of v n 
and £„. Given the nonlinear nature of the hydrodynamic 
equations of motion, a linear relation between v n and e n 
is a priori not obvious. However, a strong linear correla¬ 
tion is seen in Fig. 0c), especially for n = 2, as in several 
other initial-state models. We find that at a given cen¬ 
trality, the ratio k n decreases as n increases, indicating 
reduced conversion efficiency for higher harmonics. As 
stated above, the shape of the overlap zone is predomi¬ 
nantly elliptic, which makes the conversion of £2 into V 2 
relatively more efficient than that for higher harmonics. 
For n > 2, the shapes are governed by small-scale struc- 
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FIG. 8: (Color online) Centrality dependence of the ratio 
a Vn /(v n ), for charged hadrons with momenta pr > 0.5GeV/c 
in the AMPT+Hydro model (solid lines) as compared with 
the ATLAS data UjJ (symbols). Also shown is the ratio 
Oe n /(£n) obtained in the AMPT model (dashed lines). 


tures which makes the conversion relatively less efficient. 
For n > 2 and peripheral collisions, the short lifetime of 
the system makes the conversion even weaker. 

Fluctuations in the initial positions of the nucleons and 
the formed partons in an event lead to an event-by-event 
distribution in the initial e n as well as in the final v n . 
Figures (3a) and (3b) show the scatter plots for u 2 ver¬ 
sus £2 and V 3 versus £3 for 20 — 25% centrality Pb+Pb 
collisions in the AMPT+Hydro model. Note that e n are 
calculated at r = r sw using Eq. ©. The scatter plots can 
be described to a good approximation by a linear relation 
(v n ) = k n e n , where fc 2 = 0.327 and £3 = 0.236. Figure 
[ 3 c) shows the probability density distribution p{v n ) for 
charged particles with pt > 0.5 GeV/c at the same cen¬ 
trality in this model and compares it with the ATLAS 
data 4l|. The shapes and magnitudes of the r > 2 and 
V 3 distributions are remarkably well reproduced by the 
model. We also show in Fig. ©c) the scaled probability 
density distribution p(£ n )/k n for the AMPT initial colli¬ 
sion geometry. These too are in excellent agreement with 
the corresponding p(v n ) distributions, except at large t> 2 , 
where the scaled distribution for n = 2 drops faster. This 
is attributable to slightly enhanced conversion coefficient 
k 2 = (^ 2 } /£2 at large £ 2 , as can be seen in Fig. ©a). This 
deviation from the linear behavior has been also observed 
within recent hydrodynamic model simulations |42;]. 

These distributions can be used to calculate the respec¬ 
tive mean and variance. In Fig. ©we present the central¬ 
ity dependence of cr Vn /(v„) and <r e>l /(£„) for n = 2,3 in 
the AMPT+Hydro calculations. These results are com¬ 
pared with the ATLAS data for cr Vn /(v n ) [411 . In the 
fluctuations-only scenario, these ratios are expected to be 
about 0.52 41], which is corroborated by the model cal¬ 
culations for n = 2 in the ultracentral limit and for n = 3 
at all centralities. With increasing impact parameter, the 
smaller number of participants causes the event-by-event 
fluctuations and hence the variance a en to increase grad- 
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FIG. 9: (Color online) Transverse-momentum dependence of 
(u^(pt)) 1 ^ 2 for charged hadrons in 20 — 30% and 30 — 40% 
central Au+Au collisions at sJsnn = 200 GeV in the 
AMPT+Hydro calculations with 77 /s = 0.08 (lines) as com¬ 
pared with the experimental data from PHENIX [43|] (solid 
symbols) and STAR juLIHI (open symbols). 


ually. Now, for n = 2, the (e 2 ) increases rapidly up 
to ~ 25 — 30% centrality, and thereafter the increase is 
somewhat slower; see Fig. ©b). Consequently, the ra¬ 
tio er E2 /(£ 2 ) shows a minimum at ~ 25 — 30% centrality 
and a slower growth rate at higher centralities. Similar 
behavior is seen for the ratio oy 2 / (u 2 ), which seems to 
support the linear hydrodynamic response, at least up to 
midcentral collisions. For peripheral collisions, the fail¬ 
ure of the linear hydrodynamic response (v n oc £ n ), as 
discussed above, leads to the disagreement between the 
two ratios. 

Finally, for the sake of completeness, we show in Fig. 
© the rms values of v n (pT) in the full AMPT+Hydro 
calculation in comparison with the PHENIX and STAR 
data for Au+Au collisions at the top RHIC energy. With 
the shear viscosity to entropy density ratio of p/s = 0.08, 
which corresponds to 33% reduction from the LHC value 
0.12 estimated above, the model describes the RHIC data 
well within the systematic uncertainties. A smaller value 
of p/s at RHIC compared to that at LHC was also found 
in previous studies 18] that employed various fluctuating 
initial conditions. 


CONCLUSIONS 

In summary, within a coupled AMPT transport plus 
(2+l)D viscous hydrodynamics approach, we have stud- 












ied various aspects of the anisotropic flow coefficients in 
heavy-ion collisions at RHIC and LHC. Fluctuations in 
both the number and positions of the participant nu¬ 
cleons and the formed partons and the subsequent par- 
ton dynamics, in the event-by-event AMPT model, pro¬ 
vide realistic initial conditions for the later viscous evo¬ 
lution of the QGP and hadronic matter until it decou¬ 
ples. We find that this model provides a good agree¬ 
ment with the transverse-momentum-dependent and in¬ 
tegrated flow coefficients v n (n = 2 — 6 ) at RHIC with 
a minimal rj/s = 1/47T and at LHC with a larger value 
rj/s = 1.5/47T at various collision centralities. 

Our other results are as follows. We have isolated the 
contributions of parton cascade to the generation of ini¬ 
tial eccentricities and the collective flow, and they are 
found to be significant. Insensitivity of the results to 
the precise value of the switchover time (within a win¬ 
dow) suggests that the partonic scatterings in the AMPT 
model drive the system to near equilibrium in less than 
~ 1 fm/c. For ultracentral collisions, the model is able 
to reproduce qualitatively the unexpected hierarchy of 
v u {pt ) observed by CMS. The model agrees with the ob¬ 
served probability distributions p(v 2 ) and p(v 3 ) at mid¬ 
central collisions. There is a hint of a nonlinear hydrody¬ 
namic response even for the harmonics n = 2, 3 at large 
eccentricities. 
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